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ABSTRACT 

We present CCD UBVRI photometry of the field of the open cluster NGC 6866. 
Structural parameters of the cluster are determined utilizing the stellar density profile 
of the stars in the field. We calculate the probabilities of the stars being a physical 
member of the cluster using their astrometric data and perform further analyses using 
only the most probable members. The reddening and metallicity of the cluster were 
determined by independent methods. The LAMOST spectra and the ultraviolet excess 
of the F and G type main-sequence stars in the cluster indicate that the metallicity 
of the cluster is about the solar value. We estimated the reddening E(B — V) = 
0.074 ± 0.050 mag using the U — B vs B — V two-colour diagram. The distance 
modula, the distance and the age of NGC 6866 were derived as ft = 10.60 =b 0.10 
mag, d = 1189 ±75 pc and t = 813 ±50 Myr, respectively, by fitting colour-magnitude 
diagrams of the cluster with the PARSEC isochrones. The Galactic orbit of NGC 6866 
indicates that the cluster is orbiting in a slightly eccentric orbit with e = 0.12. The 
mass function slope x = 1.35 ± 0.08 was derived by using the most probable members 
of the cluster. 

Key words: Galaxy: open cluster and associations: individual: NGC 6866 - stars: 
Hertzsprung Russell (HR) diagram 


1 INTRODUCTION 

Distributed along the Galactic disc, open clusters are mem¬ 
bers of the thin-disc component of the Milky Way Galaxy. 
Since the stars in an open cluster are almost at the same 
metallicity, distance and age, but from different spectral 
classes, these systems provide valuable data to study the 
Galactic structure, chemical compo sitio n, stell ar structure 
and star formation processes (|Friellll995l . l2013h . Thus, the 
structural and astrophysical parameters of open clusters 
such as the reddening, metallicity, distance and age have 
to be determined through precise photometric and spectro- 
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scopic observations. Furthermore, the stellar structure and 
evolution models exhibit considerable discrepancies at dif¬ 
ferent parts of the observational main-sequence due to adop¬ 
tion of different input physics ([Thompson et al.ll2014h . Since 
open clusters allow for direct comparisons of the predictions 
of theoretical models with observations, they provide very 
important tools to better understand the physical reasons 
behind some of these discrepancies. 

The most common method used to infer the astro- 
physical parameters of the Galactic and extra-galactic open 
clusters is based on fitting stellar isochrones to the ob¬ 
served colour-magnitude diagrams (CMDs) and two-colour 
diagrams (TCDs) of the clusters. In this method, astrophys¬ 
ical parameters are estimated simultaneously via Bayesian 
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statistics. However, such simultaneous statistical solutions 
based on comparison of stellar isochrones with the photo¬ 
metric observations suffer from degeneracies between param¬ 
eters, causing lar ge uncertainties in the measured redden¬ 
ing, metallicity ([King et ahl 120051 : Ide Meulenaer et all 120131 : 
I Janes et al. issir and therefore the age values. 

Different approaches were suggested to break the 
age-reddening degeneracy in the simultaneous solutions. 
The main idea behind most of these suggestions is to 
use the maximum available wavelength range, prefer¬ 
ably including at l e ast one near-infrared (NIR) band (cf . 
Anders BB e^L_ 2004- Bridzius et al. 2008: iBilir et al.1 l20ld : 


de Meulenaer et al l2013h . lAnders et, alJltooT suggest that, 
if NIR data cannot be obtained and only observations in 
four passbands are available, one of the best photometric 
band combination for a reliable parameter determination is 
UBVI. Besides these suggestions, independent (traditional) 
and reliable methods that were developed for the determina¬ 
tion of the reddening and metallicity, such as those recently 
used bv lYontan et al.l (|2015l b can be used to constrain these 
parameters. 

Here, we focus on the CCD UB VRI observations of open 
cluster NGC 6866 (l = 79°.560, b = +6°.839) since it is o ne 
of th e clusters located in the Kepler field ([Borucki et al.l 
IB) and atmospheric model parameters and radial ve¬ 
locities for a considerable number of stars in the cluster’s 
field are present in the Large Sky Area Mu lti-Object Fiber 
Spectroscopic Telescope (LAMOST: lLuo et al IB) survey. 
The availability of photometric and spectroscopic data al¬ 
low us to calculate the cluster’s precise astrophysical and 
kinematical parameters. Following independent methods we 
find reddening and metallicity of the cluster from the stars 
with high membership probabilities and then keeping these 
two constant, we derive its distance modulus and age by fit¬ 
ting stellar isochrones to the observed CMDs. Although it 
has recently been investigated photo metrically, NGC 6 866 
is a relatively less well-known cluster (| Janes et al ](B) as 
compared to the other open clusters located in the Kepler 
field (NGC 6811, NGC 6819 and NGC 6791). The colour 
excesses, distance moduli, distances and ages obtained for 
NGC 6866 in the previous studies are summarized in Table 
1. Almost all of thes e studies have assumed a solar metal¬ 
licity for the cluster. iMermilliod et al.1 (200 8) measured the 
radial velocities of two giant stars located in the cluster field 
and estim ated a mean radial velocity o f V r = +13.68 ± 0.09 
km s" 1 . iFrinchabov &; Maiewskil (2008) derived a radial ve¬ 
locity of V r = +12.18 ± 0.75 km s -1 , a value i n agree¬ 
ment with the mean velocity given by IMermilliod et all 
( 2008 b They also found the proper motion components as 
Ha cos S = 5 .52 + 1. 17 and ns — 7 .97 + 1.09 mas yr -1 for 
NGC 6866. lMolenda-Zakowicz et~aT. ( 2009 ) detected 19 vari¬ 
able stars in the cluster field, changing V apparent mag¬ 
nitudes from 12.1 to 17.4 and B — V colours from 0.3 to 
2.2 mag. ljoshi et al.l ( 2012 4 found 28 periodic variables with 
V apparent magnitudes between 11.5 and 19.3 and B — V 
colours between 0.3 and 2.1 mag. Nineteen of these systems 
were newly identified. 

In this work, our main goal is to overcome in part 
the parameter degeneracy in the statistical solutions of the 
astrophysical parameters of NGC 6866. In Section 2, we 
summarize the observations and reductions. We present the 
CMDs, structural parameters of NGC 6866, and the mem- 


Table 1. Colour excesses ( E(B — V )), distance moduli (/x), dis¬ 
tances (d) and ages (t) collected from the literature for the open 
cluster NGC 6866. References are given in the last column. 


E{B - V) 
(mag) 

(mag) 

d 

(pc) 

t 

(Myr) 

Refs. 

0.14 

10.82 

1200 

— 

1, 2 



- 

230 

3 

0.10 

11.10 


- 

4 




650 

5 

0.17 

1 1,7.7. 

1 100 

480 

6 

0.14 

11.13 

1380 

560 

7 

0.10 

11.15 

1470 

630 

8 

0.19 

11.08 

1650 

800 

9 

0.16 

10.98 

1250 

705 

10 
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(I2014T) . 


bership probabilities of the stars in the cluster field in Sec¬ 
tion 3. In section 4, we measure the astrophysical parameters 
of the cluster. Section 5 discusses the results of the study. 


2 OBSERVATIONS 
2.1 Photometry 

CCD UBVRI images of th e ope n cl uster NGC 6866 and 
standard stars selected from lLandoltl ( 2009 4 were acquired 
on 18th July 2012 using a lm Ritchey-Chretien telescope 
(T100) located at the TUBiTAK National Observatory 
(TUG)0 in Bakirlitepe, Antalya/Turkey. With a relatively 
large CCD camera which has a field of view of about 
21.5x21.5 arcmin, this telescope is very useful for open clus¬ 
ter observations. Since the observations of NGC 6866 (Fig. 
1) and NGC 6811 were performed in the same night, the de¬ 
tails of the obs ervations and photometric reductions can be 
found in lYontan et al.l ( 2015 1. where the photometric anal¬ 
yses of NGC 6811 are discussed. Below, a brief summary is 
given. 

Short and long exposure images were obtained in each 
filter of the cluster’s field in order to cover the widest possi¬ 
ble flux range. The night was moderately photometric with 
a mean seeing of l /7 .5. A log of observations is given in Ta¬ 
ble 2. IRAIo routines were used for pre-reduction processes 
of images and transforming the pixel coordinates of the ob¬ 
jects identified in frames to equatorial coordinates. The in¬ 
strumental magnitudes of the standard stars were measured 
utilizing the aperture photometry packages of IRAF. Atmo¬ 
spheric extinction and transformation coefficients for the ob¬ 
serving system were d etermined from the observations of the 
st andard stars (lLandoltl l2009h through the equations given 
in I Janes &; Hoal (|2013l 4. which are a set of calibration equa¬ 
tions to transform instrumental magnitudes to the Johnson- 
Cousins magnitude system. For the list of the coefficients for 

1 www.tug.tubitak.gov.tr 

2 IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories. 
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Figure 1. A 30 second R -band image of NGC 6866 obtained with 
T100 telescope of the TUBITAK National Observatory. The field 
of view is about 21x21 arcmin (North top and East left). 


Table 2. Log of observations, with exposure times for each pass- 
band. N denotes the number of exposure. 


Filter 

Long Exp. time (s)xN 

Short Exp. time (s)xN 

U 

180x3 

90x3 

B 

90x3 

20x3 

V 

60x3 

5x3 

R 

30x3 

4x3 

I 

30x3 

4x3 


that particular night, see Table 2 of lYontan et al.l (|2015l ). 
Source Extractor (SExtr actor and isophotal photometry 
(jBertin fc Arnoutsl 119961 ) were performed for the cluster’s 
field. Transformation of instrumental magnitu des to stan¬ 
dard magnitudes was done as described in lYontan et al.l 

feoish . 


2.2 Spectroscopy 

Spectra of the stars in the direction of N GC 686 6 were ob¬ 
tained as part o f the LAMOST sur v ey dLuo et a D jo ij 
Cui et all [20121 : iDe Cat et al.l I20l4 l2015al bl: iRen et ahl 
20151 ). The LAMOST, also called the Guo Jing Telescope, 
is located at the Xinglong observatory, China. The tele¬ 
scope has an aperture of 4m. Its focal plane, which has a 
5° field of view in diameter, is covered with 4000 optical 
fibers connected to 16 sets of multi-obje ctive optical spec¬ 
trometers with 250 optical fibers each (IWang et a l. Il996l : 
IXing et ahlll998l ). The spectrometer obtains low resolution 
spectra (R ~ 1800) in two wavelength regions, 3700-5900 
and 5700-9000 A, using two CCD cameras. LAMOST can 
take spectra of all objects in the field of view down to a mag¬ 
nitude of V ~ 17.8 simultaneously. Its capability of tracking 
the motion of celestial objects during about four hours while 
they are passing the meridian makes the LAMOST an ideal 
instrument to collect low-resolution spectra of the objects 
in the Kepler field. Spectra obtained with the LAMOST 

3 SExtractor: Software for source extraction. 


are wavelength calibrated and intensity norm alized follow¬ 
ing the procedures detailed in I Luo et~al . (2012). For detailed 
information abou t the LAMOST, see ICui et ahl ( 2012 ) and 
Zhao et a l ] (I2012I) . 

We scanned the LAMOST database in order to find the 
stars with spectra, which are located in our field of view, and 
found 31 stars for which atmospheric parameters and radial 
velocities were measured. We presented the KIC names and 
equatorial coordinates of these stars in Table 3 (see also, 
Fig. 7). 

The stellar atmospheric parameters, e.g. effective tem¬ 
perature T e ff, surface gravity log#, and metallicity [Fe/H], 
are estimated via the stellar parameters pipeline of the 
LAMOST. The spectral data reduction procedure consists 
of three stages: i) LAMOST 2D pipeline works on reduc¬ 
tion of CCD data including dark and bias subtraction, flat 
field correction, spectra extraction, sky subtraction, wave¬ 
length calibration, sub-exposure merging and wavelength 
band combination, ii) LAMOST ID pipeline identifies spec¬ 
tral type classes and measures radial velocities for star or 
the redshifts for galaxy or QSO using cross correlation. 
The pipeline produces four prime classifications, namely 
STAR (with released subclass), GALAXY, QSO and UN- 
KNOWN, m) LAMOST Stellar Parameter pipeline (LASP; 
IWu et all 120151 : ILuo et aD 12015). automatically derives the 
stellar parameters for late A and FGK type stellar observa¬ 
tion from LAMOST survey stars. The LASP con secutively 
adopts correlation function interpolation (CFI; |Pu et ahl 


201 2P method and U niversite de Lyon Spectroscopic analy¬ 


sis So ftware fULvSS: lKoleva et al.ll2009l : IWull2009l : IWu et ahl 
l201lD method to determine the stellar parameters. The 
LASP first uses CFI method to get an initial coarse guess 
for T e //, logg, [Fe/H] that serve as the starting values for 
ULySS, then ULySS determines the parameters by minimiz¬ 
ing the x 2 values between the observation and the template 
spectra which were generated by an interpolato r b uilt o n 
the e mpirical ELODIE stellar library (|Prugniel fc SoubiranJ 
l200lD . Details for th e spectr al analy sis and propagated er¬ 
rors can be found in lLuo et al ] (I 2015 I) . 


3 DATA ANALYSIS 

3.1 Identification of stars and photometric errors 

We identified 2096 sources in the field of NGC 6866 and con¬ 
structed a photometric and astrometric c atalogue . The stel¬ 
lar ity index (SI) provided by SExtractor (Ber tin &; Arnoutsl 
Il996h was used to detect non-stellar objects, most likely 
galaxies, in our catalogue. The SI has va lues between 0 and 

1 . According to iBertin fe Arnoutsl ([ 19961 ) , a source with the 
SI close to 1 is a point source (most probably a star), while 
an extended object has an SI close to zero. A pplication of 
the SI to t he stellar fields can be found in lAndreuzzi et al.1 
1I2002I) and lKaraali et all |2004) who showed that the objects 
with an SI smaller than 0.8 can be assumed to be extended 
objects. Thus we adopted that the objects with an SI larger 
than 0.8 are most probably stars. The V apparent magni¬ 
tude versus SI diagram of 2096 objects is plotted in Fig. 

2. Resulting catalogue contains 2089 stars. Individual stars 
in the final photometric catalogue are tabulated in Table 4. 
The columns of the table are organized as ID, equatorial 
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Table 3. ID numbers, Kepler Input Catalogue (KIC) names and equatorial coordinates of the stars with LAMOST spectra in the 
direction of NGC 6866. Apparent magnitudes (V), colour indices (B — V), observing dates and signal-to-noise ratio (SNR) of the spectra 
are also given. 


ID 

KIC 

O2000.0 

(hh:mm:ss.ss) 

<52000.0 

(O , 

V 

(mag) 

B-V 

(mag) 

Date 

SNR 

1 

8329578 

20:03:00.39 

+44:16:54.13 

13.519 

0.728 

25.10.2013 

78 

2 

8329629 

20:03:04.09 

+44:14:15.65 

12.588 

0.267 

25.09.2013 

186 

3 

8128952 

20:03:04.83 

+43:59:58.60 

13.786 

0.568 

25.09.2013 

74 

4 

8263825 

20:03:12.47 

+44:06:59.03 

13.445 

0.842 

25.10.2013 

87 

5 

8329888 

20:03:22.22 

+44:16:57.95 

12.162 

0.295 

25.10.2013 

182 

6 

8329894 

20:03:22.81 

+44:15:50.46 

11.578 

1.030 

25.09.2013 

166 

7 

8264037 

20:03:26.11 

+44:10:05.45 

11.755 

0.356 

25.09.2013 

208 

8 

8264075 

20:03:28.33 

+44:07:55.26 

13.784 

0.393 

25.10.2013 

66 

9 

8264148 

20:03:32.47 

+44:06:46.17 

14.032 

0.479 

25.09.2013 

60 

10 

8330092 

20:03:34.92 

+44:14:50.16 

13.548 

0.439 

25.09.2013 

89 

11 

8197368 

20:03:44.34 

+44:05:39.30 

12.777 

0.475 

25.10.2013 

137 

12 

8330251 

20:03:45.04 

+44:15:19.86 

13.321 

0.688 

25.09.2013 

102 

13 

8129588 

20:03:45.24 

+43:59:22.51 

13.207 

0.341 

25.10.2013 

114 

14 

8197440 

20:03:48.06 

+44:02:48.38 

12.943 

0.674 

25.10.2013 

137 

15 

8264534 

20:03:54.16 

+44:06:46.04 

12.628 

0.222 

25.09.2013 

169 

16 

8264581 

20:03:57.12 

+44:08:16.86 

13.509 

0.436 

25.09.2013 

103 

17 

8330453 

20:03:57.80 

+44:16:32.51 

12.804 

0.738 

25.09.2013 

131 

18 

8264617 

20:03:59.33 

+44:10:25.84 

13.975 

0.439 

25.09.2013 

79 

19 

8264674 

20:04:02.84 

+44:11:55.63 

11.164 

0.251 

25.09.2013 

328 

20 

8330543 

20:04:03.64 

+44:15:34.07 

13.219 

0.551 

25.10.2013 

99 

21 

8264698 

20:04:03.95 

+44:10:20.72 

12.356 

0.318 

25.09.2013 

197 

22 

8197761 

20:04:09.30 

+44:04:16.39 

10.639 

0.338 

25.09.2013 

329 

23 

8396247 

20:04:09.36 

+44:19:10.03 

13.835 

1.852 

25.09.2013 

36 

24 

8330778 

20:04:16.17 

+44:12:04.66 

13.468 

0.379 

25.10.2013 

118 

25 

8330790 

20:04:16.73 

+44:12:43.66 

13.143 

0.314 

25.09.2013 

161 

26 

8264949 

20:04:18.36 

+44:09:52.47 

12.122 

0.302 

25.09.2013 

235 

27 

8265068 

20:04:25.50 

+44:10:16.54 

12.029 

0.333 

25.09.2013 

282 

28 

8198114 

20:04:32.20 

+44:05:12.84 

13.728 

0.660 

25.09.2013 

71 

20 

8265356 

20:04:43.26 

+44:08:02.05 

12.367 

0.611 

17.10.2013 

93 

30 

8331290 

20:04:44.09 

+44:15:05.77 

14.059 

0.691 

25.09.2013 

54 

31 

8265377 

20:04:44.37 

+44:08:58.53 

13.170 

0.315 

25.09.2013 

117 


Table 4. Photometric and astrometric catalogue for the open cluster NGC 6866. The complete table can be obtained electronically. 


ID <32000 J 2000 V 

(hh:mm:ss.ss) (dd:mm:ss.ss) (mag) 


U — B B — V V — R R — I /la cos 5 /is P 
(mag) (mag) (mag) (mag) (mas yr -1 ) (mas yr -1 ) (%) 


1 

2 

3 

4 

5 


20:02:57.52 

20:02:57.53 

20:02:57.58 

20:02:57.74 

20:02:57.77 


+44:11:17.39 

+44:04:17.23 

+44:14:51.51 

+44:09:23.85 

+44:12:06.25 


18.454+0.017 

18.546+0.018 

18.319+0.015 

18.873+0.021 

16.612+0.005 


1.445+0.270 

1.163+0.182 

0.703+0.097 

0.849+0.154 

0.607+0.029 


1.381+0.036 

1.149+0.033 

1.033+0.026 

1.082+0.037 

1.002+0.008 


0.902+0.022 

0.685+0.024 

0.637+0.021 

0.623+0.029 

0.607+0.007 


0.676+0.019 

0.697+0.023 

0.643+0.021 

0.743+0.028 

0.611+0.007 


9.0+6.7 


1.2+6.7 32 


coordinates, V apparent magnitude, colours, proper motion 
components and the probability of membership. The proper 
motions of the stars wer e taken from the astrometric cata¬ 
logue of lZacharias et ahl (UCAC4: l2013h . 

The errors of the measurements in the V band and 
U — B, B — V, V — R, V — I and R — I colours are shown 
in Fig. 3 as a function of the apparent V magnitude. Mean 
errors in the selected magnitude ranges are listed in Table 
5. The errors are relatively small for stars with V < 18 mag, 
while they increase exponentially towards fainter magni¬ 
tudes. Thus, we decided to use stars with V < 18 mag for 
further analysis. With this selection, the number of remain¬ 


ing stars for analysis is 1301 in the field of NGC 6866. As ex¬ 
pected, the largest errors for a given V magnitude occurred 
in the U — B colours of the stars. 

We compared our p hoto metric mea s ureme nts with 
those of Ijoshi et al.l (|2012h and I Janes et al.l ([2014) who uti¬ 
lized the same photometric bands, i.e. UBVRI , in their 
investigations. A cross-match of our ca talogue with these 
two catalogues r esulted in 800 star s for Ijoshi et al.l (|2012h 
and 429 stars for I Janes et al ] d20l4) . Comparisons between 
our and their datasets are shown in Fig. 4 and Fig. 5. In 
these figures, values on the abscissae refer to our measure¬ 
ments, while the magnitude or colour differences in the ordi- 
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Figure 2. Stellarity index (SI) for 2096 sources observed in the 
direction to the cluster NGC 6866, as a function of V appar¬ 
ent magnitude. Dashed line represents the SI of 0.8. Only seven 
sources with SI < 0.8 are classified as non-stellar sources. 


Table 5. Mean errors of the photometric measurements for the 
stars in the direction of NGC 6866. N indicates the number of 
stars within the V apparent magnitude range given in the first 
column. 


Mag. Range 

N 

av 

au-B 

&B — V 

av-R 

1 

ES 

b 

10 < V ^ 12 

18 

0.001 

0.002 

0.002 

0.002 

0.002 

12 < V ^ 13 

39 

0.001 

0.002 

0.002 

0.002 

0.002 

13 < V ^ 14 

65 

0.001 

0.004 

0.001 

0.002 

0.002 

14 < V ^ 15 

109 

0.002 

0.011 

0.003 

0.002 

0.003 

15 < V ^ 16 

192 

0.003 

0.022 

0.005 

0.004 

0.004 

16 < V < 17 

320 

0.005 

0.039 

0.009 

0.007 

0.007 

17 < V < 18 

560 

0.009 

0.083 

0.016 

0.013 

0.012 

18 < V < 19 

608 

0.016 

0.150 

0.030 

0.023 

0.021 

19 < V ^ 20 

178 

0.027 

0.183 

0.050 

0.037 

0.035 



Figure 3. V apparent magnitude versus colour and magnitude 
errors for the stars identified in the field of the open cluster 
NGC 6866. 


nates present the differences between the our and their cat¬ 
alogues. The means and the standard deviations of the mag¬ 
nitudes and colours obtained from the differences are small 
except for the standard deviations of A u and A u-b, i- e. 
au — 0.141 a nd gjj-b — 0.133 mag in the comparison with 
the Ijoshi et~aD ([20121 ) and au = 0.108 and &u-b = 0.095 
mag in the comparison with the Jan es et al ] (|2014l ). Note 
that the large mean differences and standard deviations cal¬ 
culated for the U magnitudes and U — B colours are due to 
faint field stars beyond U ~ 17 mag, for which the data are 
not used in the derivation of t he cluster p ara mete rs. Agree¬ 
ment between our study and I Janes et al.l (l2 0 14h is much 
better than that with the study of ljoshi et al.l |2012h . 

3.2 Cluster radius and radial stellar surface 
density 

We estimated the stellar density profile of the open clus¬ 
ter NGC 6866 using stars with V < 18 mag in the field. 
The central coordinates of the cluster were assumed to 
be as given in WEBD^0 database (0:2000.0 = 20^03 m 55 s , 
£2000.0 = +44°09 / 30 // ). We then calculated the stellar den¬ 
sity in an area defined by a circle with a radius of 1 arcmin 
and centered on these coordinates. From this central circle, 
we calculated the variation of stellar density using annuli 

4 http://www.univie.ac.at/webda/ 


with width of 2 arcmin. The last annulus had a width of 3 
arcmin because of a significant decrease in the number of 
stars. The resulting stellar density profile of the cluster is 
plotted in Fig. 6. 

We fitted this density profile with the iKin 3 (Il962h 
model defined as: 

P(r) = fb9 + 1 + (r/r e y' (1) 

where r is the radius of the cluster centered at the celestial 
coordinates given above. fb g , /o and r c denote the back¬ 
ground stellar density, the central stellar density and the 
core radius of the cluster, respectively. We fitted the King 
model to the observed radial density profile for NGC 6866 
and used a y 2 minimization technique to determine fb g , /o 
and r c . The best fit to the density profile is shown with a 
solid line in Fig. 6. We estimated the central stellar density 
and core radius of the cluster, and the background stellar 
density as fo = 2.28 d= 0.02 stars arcmin -2 , r c — 3.24 d= 
0.04 arcmin and fb g = 5.33 d= 0.01 stars arcmin -2 , respec¬ 
tively_ 

Ijoshi et al.l ([2012 ) estimated the central stellar density 
and core radius of the cluster, and the background stel¬ 
lar density as fo = 5.7 d= 0.7 stars arcmin -2 , r c = 2.0 
db 0.5 arcmin and fb Q = 7.85 stars arcmin -2 , respectively, 
while Ijanes et al.l ( 20141 ) found fo — 1.65 stars arcmin -2 , 
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Figure 4. Comparison of the magnitu des and col o urs fo r the 
stars observed both in this study and in ljoshi et al . (2012 b The 
means and standard deviations of the differences are shown in 
panels. 


r c = 4.70 arcmin and fb g — 0.63 stars arcmin -2 , respec¬ 
tively. Structural parame ters in this s tudy are generally in 
agreement with those in Ijoshi et al. (boil) w herea s they 
are not in agreement with those in Janes et al.l ( 2014lb This 
is mos t likely caused by the selection made by I Janes et al.l 
<l2014h since they considered only 75 stars of 11-16 mag 
within a circular field of 3 arcmin radius from the cluster’s 
centre. 


Figure 5. Same as 


Fig.4 but for I Janes et al. ( 2014 b 



3.3 CMDs and membership probabilities 

In order to derive the parameters of NGC 6866, we used its 
CMDs in Fig. 7, i.e. V vs U-B, V vs B-V, V vs V-R and 
V vs R — I. An inspection by eye suggests that the cluster is 
rather sparse. Since most of the stars brighter than V = 17 
mag are lying along a sequence, we conclude that these stars 
are probably representing the main-sequence of the cluster. 
On the contrary, most of the stars fainter than V — 17 mag 
are field stars. In order to confirm this conclusion we selected 
a circular area of the core radius and an annulus between 
10.5 and 11 arcmin from the centre of the cluster. We found 
237 and 67 stars in the circular area and the annulus, respec¬ 
tively, and presented their positions in the V vs B — V CMD 
in Fig. 8 demonstrating that most of the stars close to the 
cluster centre are generally brighter than V «17 mag and 


Figure 6. Stellar density profile of NGC 6866. Errors were de¬ 
termined from sampling statistics: 1/y/N, where N is the number 
of stars used in the density estimation. 


lying along the cluster’s main-sequence, as indicated above. 
However, the stars in the annulus are mostly fainter than 
V ~ 17 mag and they are probably field stars. 

For relatively old clusters like NGC 6866, the red 
clump (RC) stars in their CMDs can be very useful 
to determine their distances and ages since they have 
been utilized as standard candles for distance estimates 


(i.e. 

Paczvnski & Stanek 19981; Cabrera-Lavers et al. 2005. 

2007 

: Bilir et al. 2013a; Karaali et al.l 20131:lYaz Gokce et al. 

2013 

). The RC stars in V vs B — V CMD can be identified 










































































Open cluster NGC 6866 7 



Figure 7. The CMDs for the cluster NGC 6866. Black dots represent the stars with LAMOST spectra. 


in the colour range 0.7 ^ (B — R)o ^1.2 mag and the abso - 
lute magnitude range 0 ^ My ^ 2 mag ([Bilir et al . 2013b). 
Here, (B — V) o denotes the de-reddened B — V colour. In¬ 
deed, there are stars in these colour and apparent magnitude 
ranges in Fig. 7. If they are members of the cluster, these 
RC stars can be used to confirm the distance estimation of 
NGC 6866. On the other hand, the position of another small 
group of bright and blue stars in the CMDs suggests that 
the turn-off point of the cluster lie within 11 < V < 12 mag 
and 0.2 ^ B — V < 0.4 mag. 

Although identifying some stars near the RC region and 
turn-off point of the CMDs is very promising for the analy¬ 
sis, we can not confidently use these stars for further analysis 
without knowing if they are physical members of the cluster. 
In addition, determining the likely members of the cluster 
can be very useful to clearly demonstrate the cluster’s main- 
sequence. Thus, the probabilities of the stars in the field be¬ 
ing physical members of the cluster must be estimated. For 
the membership probabilities of the individ ua l stars, we em¬ 
ployed the method given by iBalaguer-Nunez et al .1 (|l998h . 
This method takes both the errors of the mean cluster and 
the stellar proper motions into account. The basic idea of 
the non-parametric method for the cluster-field separation 
in the 2-dimensional proper motion space is the empirical de¬ 
termination of the cluster and field stars’ distributions with¬ 
out any assumption about their shape. We used the kernel 
estimation technique (with a circular Gaussian kernel func¬ 
tion) to derive the data distributions. The proper motions of 
the stars were taken from the UCAC4 catalogue. Addition¬ 
ally, we com pared the results wi t h tho se of the algorithm 
published by Llavakhishvili et all (|2006h , yielding excellent 
agreement. For this, we considered rectangular coordinates 
of the stars in the field, measured in two epochs, first of our 
observations and second the UCAC4 ones. The histogram 
of the differences efficiently discriminate the members from 
non-members. 

In order to determine the most likely members of the 
cluster we applied the following procedure. First, we selected 



Figure 8. The location of stars in the V vs B — V CMD selected 
from the two different regions in the field of NGC 6866. Gray dots 
denote the stars located in the region with r c ^ 3.2 arcmin, while 
black triangles indicate the stars within an annulus between 10.5 
and 11 arcmin from the centre of the cluster. 


only the stars in a circle with a radius of 6 arcmin whose 
centre coincides with the cluster’s centre. The selected radius 
of 6 arcmin corresponds to about two times the core radius 
of the cluster. As can be seen from Fig. 6, the field stars 
are dominant beyond this radius. The median value of the 
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Figure 9. The histogram of the membership probabilities esti¬ 
mated for the stars which are located in a circle with a radius of 
6 arcmin from the centre of NGC 6866. 


membership probabilities of these stars is found as P — 50%. 
The histogram of the probabilities is shown in Fig. 9. Second, 
in order to identify the main-sequence stars of the cluster, 
we fit ted the zero age main-sequence (ZAMS) of lSung et al l 
(|2013h for solar metallicity to the V vs B — V CMD of NGC 
6866 for 12.75 ^ V ^ 17 mag using only the stars with 
P > 50%. The faint apparent magnitude limit was assumed 
to be V = 17 mag as the proper motions could be obtained 
only for the stars brighter than this magnitude. By shifting 
the fitted main-sequence to brighter V magnitudes by 0.75 
mag, a band like region in F vs B — V CMD (see Fig. 10) 
was obtained to cover the binary stars, as well. However, we 
found visually out that the stars brighter than V = 12.75 
mag have already left the ZAMS. Thus, we conclude that 
this magnitude roughly corresponds to the turn-off point 
of the cluster. Hence, we assumed that all stars brighter 
than V = 12.75 mag, which are located in the circle defined 
above and having a probability of membership larger than 
P = 50%, are the most likely members of the cluster. Finally, 
we assumed that all stars with a membership probability 
P > 50% and located within the band-like region defined 
above are the most likely members of NGC 6866 on the 
main-sequence. With this procedure, we identified 64 stars 
which are used for further analyses. These stars are indicated 
with black dots in Fig. 10. 


4 DETERMINATION OF THE 

ASTROPHYSICAL PARAMETERS OF 
NGC 6866 

4.1 The reddening 

Before the determination of the metallicity of the cluster 
using photometric observations, its reddening should be es¬ 
timated. Thus, for the determination of the colour excesses 
E(U — B) and E(B — V), we used the 64 probable members 
of the cluster selected according to the procedure in Section 
3.3. We compared the positions of these stars i n the U — B 
vs B — V TCD with the ZAMS of ISung et al.l (l2013h with 
a solar metallicity. In order to d o this, we shifted the de- 
reddened main-sequence curves of ISung et al.l ( 2013) within 
the range 0 ^ E(B — V) ^ 0.20 mag with steps of 0.001 mag 
until it fits well with the U — B vs B — V TCD of NGC 6866. 



Figure 10. V vs B—V CMD of NGC 6866 constructed using stars 
which are located in a circle of 6 arcmin radius from the ce ntre o f 
the cluster. Solid lines represent the ZAMS of lSung et al.l ( 20130 
and the one shifted by an amount of 0.75 mag to the bright V 
magnitudes. Black dots denote the most probable cluster stars 
that are identified using a procedure explained in the text. Open 
triangles indicate the stars with LAMOST spectra. 


The shift in the U — B axis was calculated b y adopting the 
following equation (Hil tner &; Johnsonlll956h : 


E(JJ — B) — E(B -V)x [0.72 + 0.05 x E(B - V)]. (2) 

To define the goodness of the fit we adopted the mini¬ 
mum x 2 method. Fig. 11 shows the U — B vs B — V TCD 
of NGC 6866 for the 64 probable members of the cluster. In 
Fig. 11 , the d a shed red fine represents the reddened ZAMS 
of lSung et al.l (12013 ") while the black dotted fine shows the 
de-reddened ZAMS of the cluster and green fines the ±1 a 
deviations. This method gives the following colour excesses: 
E(U — B) — 0.054 ± 0.036 and E(B — V) = 0.074 ± 0.050 
mag. The errors indicate the ±lcr deviations. In order to 
evaluate the E(V — R ) and E(R — I) colour exc esses we 
usedthe following equations of ICardelli et al.l (1 98,9 ): 


E(V — R) = 0.65 x E(B - V ), (3) 

E(R-I) = 0mxE(B-V). 

The colour excesses calculated are E(V — R) = 0.048±0.050 
and E(R - I) = 0.044 zb 0.050 mag. 
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B-V 

Figure 11. U — B vs B — V TCD for the main-sequence stars with 
12 ^ V ^ 17 mag of NGC 6866. The reddened and de-reddened 
main-sequence curves (ISung et al.ll2013f) fitted to the cluster stars 
are shown with red dashed and black dotted lines, respectively. 
Green lines represent ±1 a deviations. 


4.2 Metallicity and radial velocity from the 
LAMOST spectra 

Stars with the LAMOST spectra which are located in a circle 
of 6 arcmin radius from the centre of the cluster were used 
to derive the mean radial velocity and the metallicity of the 
cluster. The effective temperatures T e //, surface gravities 
logg, metallicities [Fe/H] and radial velocities V r of these 
stars are listed in Table 6. These stars are also indicated 
with green circle in Fig. 10. 

Median metallicities and radial velocities of 14 stars in 
Table 6 are —0.09 ±0.11 dex and 9.86 ± 26.01 km s -1 , re¬ 
spectively. Six stars out of 14 have the membership prob¬ 
abilities larger than 50%. Considering these six stars, we 
conclude that the median metallicity and the radial ve¬ 
locity of NGC 6866 are [Fe/H] = —0.10 ± 0.13 dex and 
V r = 10.58 ± 31.83 km s -1 , respectively. Radial velocity 
of the cluster measured in thi s study is in agreement with 
the radial velocities derived by iMermilliod et al.1 ([2008) and 
iFrinchabov &; Maiewskil (|2008l h where using only two stars 
in each study, V r = ±13.68 ± 0.09 and V r = ±12.18 ± 0.75 
km s _1 are given, respectively. We note that a comparison 
revealed that there are no common stars between our sample 
and the stars used in the other two studies. 


Table 6. Results of the spectroscopic analysis of the stars in 
the direction to NGC 6866. Membership probabilities (P), effec¬ 
tive temperatures (T e ff), surface gravities (logg), metallicities 
([Fe/H]) and radial velocities (V r ) are given. 


KIC 

p 

(%) 

T eff 

(K) 

log 9 
(cgs) 

[Fe/H] 

(dex) 

V r 

(km s _1 ) 

8197761 

00 

7153±77 

4.02±0.32 

-0.04±0.09 

-29.12±24.47 

8264534 

27 

8270±83 

3.87±0.34 

-0.03±0.10 

-6.48±20.51 

8264674 

36 

8232±59 

3.86±0.32 

-0.12±0.08 

9.47±20.87 

8264698 

61 

7730±64 

3.98±0.35 

-0.15±0.13 

12.92±32.76 

8264949 

00 

7726±53 

3.89±0.34 

-0.16±0.11 

11.59±29.85 

8264148 

40 

6977±152 

4.09±0.33 

-0.06±0.16 

14.10±22.80 

8264581 

90 

7337±113 

3.66±0.38 

0.62±0.10 

10.25±17.00 

8264037 

00 

7244±92 

3.97±0.36 

-0.04±0.11 

5.29±29.05 

8330790 

94 

7669±69 

4.05±0.33 

-0.09±0.13 

-22.79±32.92 

8264617 

76 

7177±139 

4.11±0.32 

-0.20±0.16 

10.90±25.28 

8265068 

75 

7570±50 

3.98±0.31 

-0.11±0.11 

13.59±30.90 

8264075 

17 

7177±152 

4.11±0.34 

-0.12±0.17 

10.60±26.74 

8197368 

20 

6622±107 

4.24±0.35 

-0.08±0.11 

4.79±22.63 

8330778 

53 

7268±130 

3.99±0.47 

-0.09±0.19 

8.45±59.52 


4.3 Photometric metallicity of NGC 6866 

We fo llowed the procedure described in lKaraali et al.l {2003, 

B) to measure the photometric metallicity of the open 
cluster NGC 6866. This method is mainly derived using F- 
G type main-sequence stars. Therefore, we selected 17 out 
of 64 stars of the cluster based on their colours (0.3 ^ (B — 

V) 0 ^ 0.6 mag) correspon ding to F0-G0 spectral type main- 
sequence stars (ICoxll2000h . 

We calculated the normalized ultraviolet (UV) excesses 
of the selected stars, which are defined as the differences 
between the de-reddened (U — B)q colour indices of stars 
and the ones corresponding to the members of the Hyades 
cluster with the same de-reddened (B — V )o colour index, i.e. 

5 — (U — B)o t H — (U — B)o : s- Here, the subscripts H and S 
refer to Hyades and star, respectively. Then we normalized 
the S differences to the UV-excess at (B — U)o = 0.6 mag, 
i.e. (5o.6 • 

The (U — B )o vs (B — V)o TCD and the histogram of 
the normalized <5o.6 UV excesses of the selected 17 main- 
sequence stars of NGC 6866 are shown in Fig. 12. A Gaus¬ 
sian fit to the resulting histogram allows us to calculate the 
normalized UV excess as (5o.6 = 0.030 ± 0.005 mag, where 
the uncertainty is given as the statistical uncertainty of the 
peak of the Gaussian. In order to estimate the metallicity 
([Fe/H]) of the cluster, this^mo de value wa s evaluated in 
the following equation of lKaraali et al.l ( 201 ll ): 

[Fe/H] = —14.316(1.919)(5q. 6 - 3.557(0.285)^.6 ± 0.105(0.039). (4) 

The metallicity corresponding to the mode value for the (5o.6 
distribution was calculated as [Fe/H] = —0.013±0.002 dex. 

The theoretical stellar evolutionary isochrones use the 
mass fraction Z of ah elements heavier than helium. Thus, 
we used the following relation to transform the [Fe/H] 
metallicities obtained from the LAMOST spectr a and the 
photometry to the mass fraction Z ([Mowlavi et al. Il2012h : 
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Table 7. Colour excesses, metallicities (Z), distance moduli (/i), 
distances (d) and ages (t) estimated using four CMDs. The mean 
values are given in the last line. 


CMD 



Colour Excess 

Z 


d 

t 




(mag) 


(mag) 

(Pc) 

(Myr) 

V vs U - 

B 

E{U 

- B) = 0.054 ± 0.036 

0.0154 

10.75 ± 0.11 

1271 ± 83 

850 ± 50 

V vs B - 

V 

E(B 

- V) = 0.074 ± 0.050 

0.0154 

10.53 ± 0.11 

1148 ± 75 

800 ± 50 

V vs V - 

R 

E(V 

- R) = 0.048 ± 0.050 

0.0154 

10.48 ± 0.10 

1122 ± 68 

800 ± 50 

V vs R - 

I 

E(R 

- I) = 0.044 ± 0.050 

0.0154 

10.65 ± 0.10 

1214 ± 75 

800 ± 50 




Mean 

0.0154 

10.60 ± 0.10 

1189 ± 75 

813 ± 50 


Z 0.04 + 10-t Fe / H ]' ® 

Hence, we calculated Z = 0.010 and Z = 0.012 from the 
[Fe/H] metallicities obtained from the LAMOST spectra 
and the photometry, respectively. Since these abundances 
are very close to each ot her an d to t he solar v alue, which is 
given as Z — 0.0154 bv iBressan et al ] <l2012h . we prefer to 
use the solar abundances in the determination of the astro- 
physical parameters of the cluster. 


4.4 Distance modulus and age of NGC 6866 

The reddening, metallicity, distance modulus and age of 
a cluster can be simultaneously determined by fitting the 
theoretical stellar evolutionary isochrones to the observed 
CMDs. In our case, we measured metallicity and reddening 
of the cluster using reliable traditional methods, as noted 
above. Then, in order to derive its distance modula and 
age simultaneously, we fitted the CMDs of NGC 6866 with 
the theoretical isochr ones provided by th e PARSEC syn¬ 
thetic stellar lib rary (| Bressan~et ID l2012h , which was re¬ 
cently updated (|Tang et alJ 120141 : IChen et al.ll2014h . Since 
we keep the metallicity and reddening of the cluster as con¬ 
stants during the fitting process, it is expected that the 
degeneracy/indeterminacy of the parameters will be less 
than that in the statistical solutions with four free astro- 
physical parameters. In Fig. lj^ th e best fitted theoretical 
isochrones given by IBressan et al.1 (2012b for Z — 0.0154 
and t = 800 — 850 Myr are overplotted in the CMDs. The 
estimated astrophysical parameters of NGC 6866 obtained 
from the best fits to the CMDs are given in Table 7. Errors 
of the parameters are derived by visually shifting the theo¬ 
retical isochrones to include all the main-sequence stars in 
the observed CMDs. 


4.5 Galactic orbit of the cluster 


In order to estimate the parameters of the Galactic or¬ 
bit of the cluste r, we f oll owed the procedu re described 
in iDine scu et all (|l999h , ICoskunoglu et all (|2012h and 
iBilir et alJ (|2012h . We first performed a test-particle inte¬ 
gration in a Milky Way potential which consists of a loga¬ 
rithmic halo, a Miyamoto-Nagai potential to represent the 


Galactic disc and a Hernqu ist potential to model t he bulge. 
We then used values from ICoskunoglu et al.1 i2011 ) for LSR 


corrections. 

We calculated Galactic orbits of the six cluster stars 
with a membership probability larger than 50% for which 



(B-V) 0 





Figure 12. The (U — B )o vs (B — V )o TCD (upper panel) and the 
histogram (lower panel) for the reduced UV-excess for 17 main- 
sequence stars used for the metallicity estimation of NGC 6866. 
The solid line in the upper panel represents the main-sequence 
of Hyades cluster while the one in the lower panel shows the 
Gaussian fit of the histogram. 


LAMOST data are given in Table 6. Median values of the 
radial velocity and proper motion components of these stars, 
and the mean distance of the cluster were taken as the in¬ 
put parameters for the cluster’s Galactic orbit estimation: 
V r — 10.58 km s _1 , (i a cosS — —3.30 and (is — —5.65 mas 
yr -1 , and d — 1189 pc, res pectively . T he pr oper motion com¬ 
ponents were taken from IZacharias et al • (1201 3) , while the 
radial velocities of the cluster stars and their distances were 
found in this study. Proper motions used in our study are al¬ 
most the same as used bv lWu et al. ( 2009 ) ((i a cos 5 — —3.33 
and (is — —5.03 mas yr _1 ) who also calculated the parame¬ 
ters of the Galactic orbit of the cluster. Galactic orbits of the 
stars were determined within an integration time of 3 Gyr 
in steps of 2 Myr. This integration time corresponds to min¬ 
imum 12 revolutions around the Galactic center so that the 
averaged orbital parameters can be determined reliably. In 
order to determine the Galactic orbit of the cluster, means 
of the orbital parameters found for the stars were adopted 
as the orbital parameters of the cluster. 

In Fig. 14, representations of Galactic orbits calculated 
for the cluster stars and the cluster itself are shown in the 
X — Y and X — Z planes. Here, A, Y and Z are heliocentric 
Galactic coordinates directed towards the Galactic centre, 
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U-B B-V V-R R-l 


Figure 13. Four CMDs for the stars located in a circle of 6 arcmin radius from the center of NGC 6866. The most probable members 
of the cluster are indicated with black circles. These stars are fitted to the isochrone determined in this study (solid line). The dashed 
lines indicate the isochrones with estimated age plus/minus its error. 


Galactic rotation and the north Galactic pole, respectively. 
The cluster’s apogalactic ( Rmax ) and perigalactic ( Rmin) 
distances were obtained as 9.78 and 7.76 kpc, respectively. 
The maximum vertical distance from the Galactic plane is 
calculated as Zm ax = 160 pc. When determining the ec¬ 
centricity projected on to the Galactic plane, the following 
formula was used: e — ( Rmax - Rmin)/(Rmax + Rmin)- The 
eccentricity of the orbit was calculated as e = 0.12. This 
value shows that the cluster is orbiting the Galaxy with a 
period of P or b = 156 Myr as expected for the objects in the 
solar neighborhood. 

4.6 Luminosity and mass functions of the cluster 

The luminosity function (LF) is defined as the relative num¬ 
ber of stars in the unit absolute magnitude range. Since 
the correction for the non-member stars in a cluster field is 
very important in the estimation of the luminosity function, 
we decided to demonstrate the effect of non-member stars 
using the following procedure. First, we selected the main- 
sequence stars with V ^ 17 mag located in a circular field of 
6 arcmin radius from the centre of the cluster. Note that the 
stars fainter than V ~ 17 mag h ave no accurate astr ometric 
data in the UCAC4 catalogue (|Zacharias et al.ll2013l ). Thus, 
the membership probabilities can not be estimated for the 
stars fainter than V ~ 17 mag. Then we removed the stars 
located out of the band-like region presented in Fig. 10. Ad¬ 
ditionally, we also selected the stars located near the turn-off 
point of the cluster and brighter than V = 12.75 mag. This 
selection procedure resulted in 116 stars with the member¬ 
ship probability P > 0%. For the main-sequence stars, the 
apparent magnitude V ^ 17 corresponds to a mass range of 
3.2 > M/M© > 0.8 for the cluster. In order to demonstrate 



-12 -8 -4 0 4 8 12 

X (kpc) 

Figure 14. The Galactic orbital motions (grey dashed lines) of 
the six cluster stars with P > 50%, for which LAMOST spec¬ 
tra are available, in the X — Y (a) and X — Z (b) planes. The 
cluster’s mean orbit is indicated with a blue line. The black plus, 
red circle and green triangle symbols in panel (a) represent the 
Galactic centre, and current locations of the Sun and NGC 6866, 
respectively. 
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M v 

Figure 15. The luminosity functions of NGC 6866 estimated 
for the stars with the membership probability P ^ 20% (a) and 
P ^ 50% (b). 



log M/M 0 

Figure 16. The mass functions of NGC 6866 estimated for the 
stars with the membership probability P ^ 20% (a) and P ^ 50% 
(b). 


the effect of non-member field stars, the LFs of the cluster 
estimated for the stars with the membership probabilities of 
P ^ 20% (N = 87) and P ^ 50% (N = 64) are displayed in 
Fig. 15. The LFs of NGC 6866 in Fig. 15 indicate that the 
LF of the cluster has a maximum value towards Mv ~ 5.5 
mag. 

The mass function (MF) denotes the relative number 
of stars in a unit range of mass centered on mass M. It 
represents the rate of star creation as a function of stellar 
mass. Theoretical m odels provid ed by the PARSEC syn¬ 
thetic stellar library (|Bressan et al.|[201^ i were used to con¬ 
vert the LFs to MFs for NGC 6866. Resulting MFs are pre¬ 
sented in Fig. 16. The slope x of mass function can be de¬ 
rived by using the following linear relation: log (dN/dM) = 
— (1 + x) log(M) + C, where dN represents the number 
of stars in a mass bin dM with central mass of M, and 
C is a constant. The slopes of the MFs are found to be 
x = 1.48 dz 0.21 and x =s 1.35 d= 0.08 for the stars with the 
membership probabilities of P ^ 20% and P ^ 50%, re¬ 
spectively. Since these values are in agreement in errors, we 
conclude that the effect of non-member stars can be negligi¬ 
ble. As the stars with the membership probability P ^ 50% 
constitute a more pure sample for the cluster stars, we adopt 
the MF slope x = 1.35 ± 0.08 for NGC 6866. This MF slope 
is i n a perfect agreement with the value of 1.35 given by 
ISalpeterl ( 19551 ) for the stars in the solar neighbourhood. 


5 DISCUSSION 

The reddening, metallicity, distance modulus and age of 
a cluster can be simultaneously determined by fitting 
the theoretical stellar evolutionary isochrones to the ob¬ 
served CMDs. However, astrophysical parameters derived 
from the simultaneous solutions suffer from the reddening- 


age degeneracy (cf . [Anders et al. 2004: King et all [2005; 
iBridzius et al.|[2008l : Ide Meulenaer et al.ll2013h . Thus, inde¬ 
pendent methods developed for the determination of the as¬ 
trophysical parameters are very useful to reduce the number 
of free parameters. In this study, using such methods, the 
reddening was inferred from the U — B vs B — V TCD of 
the cluster and found as E(B — V) = 0.074 d= 0.050 mag 
and the metallicities estimated from the spectroscopic and 
photometric observations were found to be roughly in agree¬ 
ment and assumed to be the solar value. Hence, keeping 
these two parameters as constants, we derived the distance 
modula and age of NGC 6866 by fitting its observed CMDs 
to the theoretical isochrones as n — 10.60 ±0.10 mag and 
t = 813 d= 50 Myr, respectively. With this method, we at¬ 
tempt to break in part the reddening-age degeneracy. 

Although the same metallicity values were adopted 
in the determination of the astrophysical parameters of 
NGC 6866, a comparison of Tables 1 and 7 reveals that 
the reddening, the distance modula, and the distance in this 
study were found to be smaller than those reported in the 
previous studie s. The determined age of the clust e r is i n 
agreement with lCiines et al] (|2012h and I Janes et al ] ll2014h . 
Si nce our photometr ic dat a is in good agreem ent with those 
of Ijoshi et al.l ( 20121 ) and I Janes et al.l (120141 ), we conclude 
that the main reason for the disagreement of the results can 
be due to the chosen isochrone and age determination differ¬ 
ence. It would be confident to conclude that our results are 
reliable since a more comprehensive approach which tries 
to minimize degeneracies among parameters is used in this 
study. 

There is only one star with P > 50% in the RC region of 
the CMDs of the cluster. The apparent magnitude and the 
colour index for this star are V — 11.699 and B — V — 0.979 
mag, respectively. We adopt E(B — U) — 0 .074 ± 0.050 (see 
Section 4.1) and Mv — 1.0 =b 0.2 mag (jBilir et aL .1 2013b ) as 


































































the intrinsic colour excess and absolute magnitude of this RC 
star, respectively. Using the Pogson’s relation, the distance 
of the RC star of the cluster is calculated as d— 1241 ±110 
pc, in agreement with the mean distance given in Table 7, 
d = 1189 ± 75 pc. 


Using the most likely members of the cluster for which 
the radial velocities are taken from the LAMOST database, 
we estimated the Galactic orbit of NGC 6866. The eccen¬ 
tricity projected on to the Galactic plane for its orbit is 
calculated to be e = 0.12, which implies that the cluster 
is located in the thin-disc component of the Galaxy. Galac¬ 
tic orbitaDparameters of the cluster were also calculated 
by Wu et al] (l2009|j usi ng the mean radial velocity given in 
iMermilliod et al.f([2008l j . They calculated the eccentricity of 
the orbit of the cluster to be e = 0.08 and the maximum 
vertical distance from the Galactic plane Z macc = 220 pc. A 
comparison shows that the eccentricity found in o ur study 
is about 50% larger than that calculated in lWu et al. (2009) 
probably since value in this study is about 30% smaller 
than their ones. 

The luminosity function of NGC 6866 has a maximum 
value near My ~ 5.5 mag. The derived slope of the mass 
function for the cluster is x = 1.35 ± 0.08 w hich is in a good 
agreement with the value of 1.35 given bv ISalpeterl (|l955l j 
for the stars in the solar neighbourhood. 
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